Closure of two dimensional turbulence: the role of pressure gradients 
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Inverse energy cascade regime of two dimensional turbulence is investigated by means of high 
resolution numerical simulations. Numerical computations of conditional averages of transverse 
pressure gradient increments are found to be compatible with a recently proposed self-consistent 
Gaussian model. An analogous low order closure model for the longitudinal pressure gradient is 
proposed and its validity is numerically examined. In this case numerical evidence for the presence 
of higher order terms in the closure is found. The fundamental role of conditional statistics between 
longitudinal and transverse components is highlighted. 
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The existence of two simultaneous inertial ranges in two- 
dimensional turbulence, as a consequence of coupled en- 
ergy and enstrophy conservation, is one of the most im- 
portant phenomena in statistical fluid mechanics jjj. At 
variance with 3-D-turbulence, the energy injected into the 
system at scale if flows toward the large scales, while the 
enstrophy cascades down on the small scales. Because of 
the inverse energy cascade, the Navier-Stokes equations, 



dtUt 



UjdjUi = 



-dip + vd Ui + ft 



(1) 



which rule the evolution of an incompressible {dim = 0) 
velocity field, cannot reach a steady-state unless an en- 
ergy sink at large scales is added. Alternatively one can 
consider an ensemble of solutions of (|l|) with a fixed en- 
ergy value below the condensation level 0, i.e. with 
an integral scale L(t) (growing in time as i 3 / 2 ) still 
much smaller than the system size. Because of the scal- 
ing of the characteristic times, the small scales (iner- 
tial range) in the system If <C r <C L can be consid- 
ered in a stationary state. One of the most challeng- 
ing problems is to understand the statistics of velocity 
fluctuations Au(r) = u(x + r) — u(x) ||. In homo- 
geneous and isotropic turbulence it amounts to study 
the joint probability density function (PDF) P(U, V, r) 
of longitudinal U and transverse V velocity differences 
where Au = USt + Vy and x = -. Recently exper- 
imental Q and numerical (H) investigations in two di- 
mensions have shown that the probability distribution 
of the pure longitudinal P(U,r) and transversal P(V,r) 
velocity differences at inertial scales display a close-to- 
Gaussian statistics with undetectable intermittency cor- 
rections to structure function exponents. Although the 
establishment of normal scaling in all inverse cascades 
seem to be generic j(| nevertheless the Gaussianity of the 
statistics in inverse cascade of the forced two dimensional 
turbulence is remained to be understood. From ([l]), a 
set of equations for generic mixed structure functions, 
i.e. S n:m (r) = (U n V m ) = A nym r^ n ' m have been obtained 
1^,1). In H those equations are elaborated from the joint 



PDF equation. Unfortunately, the PDF equation is not 
closed, resembling the well known closure problem in tur- 
bulence. In the inverse energy cascade regime, dissipative 
contributions can be neglected so the remaining unclosed 
terms are the longitudinal and transversal pressure gradi- 
ents increments. Recently Yakhot ^,|) suggested a self- 
consistent model for the pressure gradient increments and 
succeeded to obtain a Gaussian distribution for the trans- 
verse PDF, P(V,r). Although the experimental Q and 
numerical J5[ observations support the Gaussian result of 
the effective low order model, nevertheless a direct nu- 
merical computation of the pressure gradient increment 
contribution is still lacking. 

The main aim of this work is to compare the numerical 
evaluation p0[ of transverse and longitudinal components 
of pressure gradient increments with the theoretical pre- 
dictions of a recently introduced closure scheme. For the 
first time, we emphasize the importance of velocity mixed 
conditional averages such as (U\V,r) and (V 2 \U, r) gen- 
erally arising in the pure longitudinal or transversal PDF 
equations. To our surprise the existence of such objects 
has been neglected in all the previous theoretical mod- 
clings. As an essential step for the description of pure 
velocity statistics we numerically evaluate the behavior 
of these new objects for which some effective models are 
proposed. Such an investigation provides a direct check 
of the closure model. 

By standard statistical tools starting from the 

Navier-Stokes equations (|lj), it is possible to derive the 
following exact PDF equation for joint transversal and 
longitudinal velocity increments: 
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P{U,V,r) 



[e(c£ + d v ) + d v V x , u + dvTy, v ] P{U, V, r) 



(2) 



where e = (fiUi) is the rate of energy input and the 
conditional transversal, V y>v — (Ad y p\U,V,r), and lon- 
gitudinal, V x ,u = {^d x p\U, V, r), pressure gradient in- 



1 



crements are the only unclosed terms. In pure longi- 
tudinal and transversal PDF equations other unknown 
quantities play role. Indeed, by integrating (||) over 
U or V the terms (U\V,r) = UP(U\V,r)dU and 

{V 2 \U,r) = J+^V 2 P(V\U,r)dV appear in the pure 
transverse or longitudinal PDF equations respectively, 
pinpointing the statistical dependence between longitu- 
dinal and transversal components. Let us start with the 
transversal one, for which the knowledge of (Ad y p\V, r) = 
f*™ V y>v P{U\V, r)dU and of (U\V, r) is sufficient to close 
the equation. Following the recently proposed closure ||, 
we assume a second order expansion for the transverse 
pressure gradient increments V ViV in terms of local veloc- 
ity increments U and V. Even if the locality assumption 
is not based on rigorous grounds, there are some argu- 
ments in support of its plausibility ||. Once locality is 
accepted, keeping only second order terms is motivated 
from the fact that for Gaussian fields only quadratic com- 
binations of U and V appear |l2). Some physical con- 
straints simplify the expansion even further |l3|,[l4| , end- 
ing with Yakhot ansatz Q 

(Ad y p\U,V,r) = -h— -b(er) 1 / 3 -. (3) 
r r 

To directly check the closure one has to compute quan- 
tities like V y v . However, to be more quantitative, here 
we numerically compute (Ad y p\V,r) and (Ad y p\U, r) for 
which we have a better statistics. For symmetry rea- 
sons (Ad y p\U,r) = as confirmed by simulations, and 
we are left with the analysis of the term (Ad y p\V,r). 
We start by writing the quantities of interest in a scale 
invariant form. For a scale invariant solution for the 
PDF equation, i.e. P(V,r) = P(F/(er) 1 / 3 ) = P(X), 
is sufficient to require scale invariance of (U\V,r) and 
(Ad y p\V,r). We thus define {U\V,r) = {er) 1 / 3 F(X) and 
(Ad y p\V, r) = [(er) 2 / 3 /r]G(X). The major challenge now 
is to determine the functional form of G(X) and F(X). 
Taking into account the symmetries of ([!]), we assume 
for (U\V,r) an even polynomial expansion in V. Invok- 
ing the homogeneity, (U\V,r) = 0, leads to the low order 
expansion: 

F{X) =C 2 {-A Q a + X 2 ) 1 (4) 

meaning that positive (negative) longitudinal velocity in- 
crements correspond to large (small) transverse veloc- 
ity increments. Furthermore, by integrating (^|) over 
U one obtains G(X) = -hXF(X) - bX. Apparently 
this is a two parameter expansion, however the con- 
straint VP~y~^ = §,§ implies hX 2 F(X) = -bX 2 . Since 
X 2 F(X) = A\ t 2 = 1/2, one ends up with the relation 
hA\ t 2 = ~bA$ 2- The important fact is that this expan- 
sion is consistent with Gaussianity of transverse fluctua- 
tions and also gives a reasonable account for pressure con- 
tributions in the structure function equations |7|,p^5|.pii[ . 



Indeed plugging the expansion for F and G in the di- 
mensionless transverse PDF equation, one obtains the 
Gaussian result P(X) — exp(-Jf 2 /24 ,2) Hps), which is 
consistent with simulations and experiments [HJq] . Since 
positivity and finiteness of the PDF fixes the constant 
C 2 = 1/(44^2) and h = 4/3 therefore 4 2 , = 3/54 , 2 is 
the only free parameter of the theory p|Jl^]. Therefore, 
within second order approximation one has 

(Wr) = ±A ! -{-A ,2+X 2 ) , (5) 
(er) 2 / 3 fx X 3 \ 

which up to about two standard deviations agree remark- 
ably well with the numerical data (see Fig. |l|). Moreover, 
using (||) as a fitting function, we obtained 4.2,0 = 11 ± 1 
which is close, within the statistical errors, with the 
value obtained in previous experimental P,p^|, numer- 
ical P,p|]Tg|]) and analytical pl| studies. We remark that 
the good agreement of DNS data with (^|) provides a first 
evidence (even if numeric) for the plausibility of the lo- 
cality assumption. 
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FIG. 1. (a) (U\V,r) and (b) (Ad y p\V,r) computed at 
r = 0.025 (boxes) and r = 0.037 (circles). Empty symbols 
refer to the Gaussian forcing and full ones to the restricted in 
wavenumber one. The full lines represent predictions (jE|) and 
(|) With ,42,0 = 11. 
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However one can verify that assuming higher order poly- 
nomials for F(X), can result in non zero higher order 
terms in G{X). Indeed plugging the Gaussian result 
in the equation for P(X), for any order consistent with 
Gaussianity, G(X) is expressible as a functional of F(X). 
So we obtain Jill , 



G(X) = - 



X 



F(X')e 2A °.2dX' 



(7) 



Plugging the self consistent low order model (^) in (0) 
will reduce the proposal of Yakhot. Equation (Q) pro- 
vides a way to generalize the model G(X), in a self con- 
sistent way, to higher order polynomials. It is evident 
that higher order terms in F(X) can lead to higher order 
terms in G(X). So indeed one may not be able to model 
F(X) and G(X) independently, provided the Gaussian 
distribution for transverse fluctuations is assumed. 
Let us now consider longitudinal component of pressure 
gradient increment V x . u — (Ad x p\U,V,r), which has a 
major role in determining the main dynamical aspect 
of the inverse cascade, i.e. the non equilibrium energy 
flux. In contrast to the transversal case, for the longi- 
tudinal case both (Ad x p\U,r) = f+™ V x , u P(V\U,r)dV 

and (Ad x p\V, r) = V X , U P(U\V, r)dU are non-trivial. 
However, the resulting longitudinal PDF equation in- 
volves only the (Ad x p\U,r) and the velocity conditional 
average (V 2 \U,r), as one can verify by integrating @ 
over V. Therefore, only the knowledge of these two condi- 
tional averages is sufficient to close the longitudinal PDF 
equation. Again the existence of the velocity conditional 
average indicates the importance of correlation effects in 
pure longitudinal statistics. The very existence of a non- 
equilibrium flux implies that P(U,r) = P(—U, —r), hence 
the PDF equation would preserve the same invariance, 
i.e. (V 2 \ - U,-r) = {V 2 \U,r) and (Ad x p\ - U,-r) = 
— (Ad x p\U,r). Scaling invariance of the PDF equation 
implies scaling invariance of V x , u and (V 2 \U,r). Anal- 
ogous to the transversal case, we assume a local scale- 
invariant expansion for (Ad x p\U, r) and (T^ 2 |£7, r), and we 
seek for a low order closure in terms of Y = U/ (er) 1 / 3 . So 
defining (Ad x p\U,r) = [(er) 2 / 3 /r]H(Y) and (V 2 \U,r) = 
(er) 2 / 3 M(Y"), we propose the following expansion 



H(Y) = E[Y 2 - -M(Y) 



6 



-Y 



M(Y)=A 2|0 (f + ^ 



3 2A> o 



(8) 
(9) 



The coefficients of the three terms in the conditional pres- 
sure gradient are constrained by homogeneity, isotropy 
and incompressibility (i.e. YH(Y) = and H(Y) = 0). 
We observe that having reduced the expansion of M(Y) 



at the first order, the only new coefficient is the con- 
stant E. In Fig. 2 we show the numerical evaluation of 
H{Y) and M(Y). From the figure a low order expan- 
sion in terms of Y can be inferred for both these objects. 
However, concerning M(Y) the result is hardly distin- 
guishable from an almost constant value. From a best fit 
we found E = —0.39 with an error bar around 20%. If 
the longitudinal fluctuations were purely Gaussian then 
these models might be considered as a better approxi- 
mation for H(Y) and M(Y). However, the longitudinal 
statistics is just nearly Gaussian, indeed the non zero flux 
implies a non-zero skewness and to the non zero odd or- 
der structure functions S^n+i.oM — ^2n+i,o( er )^ 2n+1 ^ / ' 3 - 
Furthermore a very important observation in || in- 
dicates that the hyper-skewness of higher orders, i.e. 
S2n+i,o/S<^o + ■> increases with order and can not be 
considered as a small parameter. So the expectation from 
any kind of modeling for H(Y) and M(Y) is taking care 
of these fine details of the distribution. It seems im- 
probable to have access to these fine details within a one 
parameter low order closure or other low order models. 
As a quantitative check one can plug the low order ex- 
pansion in the longitudinal PDF equation. Then it is 
straightforward to obtain the following prediction 
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FIG. 2. (Ad x p\U,r) computed at r = 0.025 (boxes) and 
r = 0.037 (circles). Empty symbols refer to the Gaussian 
forcing and full ones to the restricted in wavenumber one. 
The full line is fitted with E — — 0.39. In the inset we show 
(V 2 \U,r) which up to two standard deviations seems to be 
constant and at larger values fluctuates. 

Substituting the numerical value of E we obtain for 
the hyper-skewness A^fi/A^ ~ 0.449 and Affi/A^Q ~ 
5.674. Comparing these numbers with the correspond- 
ing numerically obtained values, A^^/A^ q ~ 0.25 and 

A^fi/A^Q ~ 1.55, shows a large difference. The 
fourth and sixth order hyper-flatnesses calculated from 
the closure correspondingly are A^q/A^q ~ 3.29 and 
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Aj.o/Ajq ~ 20.03. Comparing to the Gaussian values 
the deviations are getting bigger with the order but still 
the errors are smaller in the even part with respect to the 
odd part of the statistics. This is an important indica- 
tion that, one has to consider higher order expansions in 
order to be consistent with higher order statistics. There- 
fore, in spite of the fairly good compatibility between the 
low order closure for H(Y) and M(Y) and their direct 
measurement in two standard deviations the fine details 
of the distribution are not recovered by them. This con- 
firms the observation in || that these fine details are 
buried in the very far tails of the antisymmetric part of 
longitudinal PDF. 

In conclusion, the dynamical role of the pressure gra- 
dient and velocity conditional averages in establishing 
the velocity increments statistics has been highlighted 
and numerically investigated. The transversal compo- 
nents of the velocity statistics has been found Gaussian, 
in agreement with previous numerical and experimen- 
tal observations. Low order expansions for the transver- 
sal conditional pressure gradient and (U\V), which have 
been proposed (in a closely related approach) in the con- 
text of a self-consistent closure have been found in 
good agreement with the DNS data up to two standard 
deviations. Furtherly, we proposed a generalization of 
the expansion which is order by order consistent with 
Gaussianity of the transverse statistics. Concerning the 
longitudinal statistics we found that the low order clo- 
sure for the conditional pressure gradient and (V 2 \U), 
although in fairly good compatibility with DNS data, 
is inconsistent with the fine details of the longitudinal 
PDF, which bear the information of the antisymmetric 
PDF tail. This indicates that unlike the transverse statis- 
tics a complete description of the longitudinal statistics 
calls for higher order terms in the expansions |16[) . It is 
worth emphasizing that these modelings are not just a 
naive fitting: the free parameters are fixed via realisabil- 
ity conditions in the PDF equations and have been tested 
numerically. Let us finally remark that the importance of 
the conditional averages goes far beyond the assessment 
of closures for two-dimensional turbulence: the impor- 
tant message is that any theoretical approach to pure 
longitudinal (transversal) velocity statistics cannot dis- 
regard the reciprocal dependence between longitudinal 
and transversal components. We consider the investiga- 
tion of such objects also in three dimensional turbulence 
a necessary step. 

We acknowledge discussions with A. Celani, 
G. Falkovich, T. Gotoh, R.J. Hill, M.R. Rahimitabar 
and V. Yakhot. G.B. has been partially supported by 
the EU under the Grant ERB FMR XCT 98-0175 "In- 
termittency in Turbulent Systems", and M.C. by the 
the EU under the Grant HPRN-CT-2000-00162) "Non 
Ideal Turbulence". J.D. acknowledges the partial sup- 
port by Deutsche Forschungsgemeinschaft (DFG). We 
also acknowledge the allocation of computer resources 



from INFM "Progetto Calcolo Parallelo" . 



[2] 

[3] 

[4] 
[5] 

[6] 

[7] 

M 



[10] 



[ii] 

[12] 
[13] 



[14] 



[15] 

[16] 
[17] 
[18] 

[19] 



R. H. Kraichnan, Physics of Fluids, 10, 1417 (1967); 
G. K. Batchelor, Phys. Fluids. Suppl. II 12, 233 (1969); 
R. H. Kraichnan, and D. Montgomery, Rep. Prog. Phys. 
43,35 (1980). 

L. Smith and V. Yakhot, Phys. Rev. Lett., 71, 352 
(1993). 

U. Frisch, Turbulence, (Cambridge University Press, GB, 
1995) 

J. Paret and P. Tabeling, Phys. Fluids 10, 3126 (1998). 
G. Boffetta, A. Celani and M. Vergassola, Phys. Rev. E, 
61, R29 (2000). 

G. Falkovich, K. Gawedzki and M. Vergassola, 

Rev.Mod.Phys. 73, 913 (2001). 

R.J. Hill, J. Fluid Mech. 434, 379 (2001). 

V. Yakhot, Phys. Rev. E, 63, 26307 (2001). 

V. Yakhot, P hys. Rev. E, 60, 5544 (1999); |arXiv:chao- 

dyn/9909017j. 



DNS of (y) have been performed by means of a fully 
dcaliascd standard pseudo-spectral code on a doubly pe- 
riodic square domain of size (2tt) 2 with 2048 2 grid points. 
Energy is injected into the system at a constant rate e by 
means of an isotropic Gaussian forcing, /, concentrated 
on small scales, if, with correlation (fi(r, t) fj (0, t')) = 
8ijS(t—t') F(r/£f) [g|. We also considered a Gaussian forc- 
ing in a restricted band of wave numbers. Viscous term in 
([j]), as customary, has been replaced by a hyper- viscous 
term (here of order eight). Time evolution is obtained 
by a second-order Runge-Kutta scheme, and integrations 
have been stopped when L(t) is still well below the sys- 
tem size to avoid condensation and averages are made 
over several independent realizations. 
A. S. Monin and A. M. Yaglom, Statistical Fluid Mechan- 
ics vol II (The MIT Press, Cambridge, Massachusetts, 
1975) 

T. Gotoh, ITP lecture notes of program on Hydrodynamic 
Turbulence, 2000. 

For P(V,r), the quantity of interest is the conditional 
average {Ad v p\V,r} = J V y , v P(U\V,r)dU. Due to rota- 
tional symmetry one has parity invariance, i.e. P{V, r) = 
P(—V,r), which implies the invariance of (Ad y p\V,r) for 
V^-V. Taking into account that (U\V, r) = (U\ -V,r), 
(Ad v p\V,r) would be an odd function of V. Moreover, 
{Ad y p\V, r)=j{Ad yP \V, r)P(V)dV=0. 
In order to guarantee the correct dimensionality, the Kol- 
mogorov scaling U~r 1//3 has been imposed in (H), i.e. we 
implicitly assumed the absence of intermittency, which 
should be consistent with the outcomes of the theory. 
S. Kurien and K.R. Sreenivasan, Phys.Rev.E 64, 056302 
(2001). 

G. Boffetta, M. Cencini and J. Davoudi, in preparation. 

J. Sommeria, J. Fluid Mech. 170, 139 (1986). 

M.E. Maltrud and G.K. Vallis, J. Fluid Mech. 228, 321 

(1991). 

P. Olla, Phys. Rev. Lett. 67, 2465 (1991). 



4 



